Stress induced stripe formation in Pd/W(110) 



OO 

o 
o 

(N 

< 

(N 



C/3 



o 



> 

O 
O 

rn 

od 

o 

00 
O 



X 



T. O. Mente§/ A. Locatelli,^ L. Aballe,^ and E. Bauer^ 

^ Sincrotrone Trieste S.C.p.A., Basovizza-Trieste 34012, Italy 
^CELLS-ALBA Synchrotron Light Facility, C3 Campus Universitat Autonoma de Barcelona, 08193 Bellaterra, Barcelona, Spain 
'^Department of Physics, Arizona State University, Tempe, Arizona 85287-1504 

(Dated: August 21, 2008) 

A stress-induced stripe phase of submonolayer Pd on W(llO) is observed by low-energy electron 
microscopy. The temperature dependence of the pattern is explained by the change both in the 
boundary free energy and elastic relaxation energy due to the increasing boundary width. The stripes 
are shown to disorder when the correlation length of the condensed phase becomes comparable to 
its period, while the condensate to lattice-gas transition takes place at a higher temperature, as 
revealed by low-energy electron diffraction. 
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It is well known that surface stress may give rise to pe- 
riodic structural modulations [H . The most common ex- 
amples are surface recontructions on clean and adsorbate- 
covered crystal faces. In the presence of competing long- 
range and short-range interactions, restructuring of the 
surface can take place at mesoscopic length scales 

Several experiments have shown the formation of meso- 
scopic patterns as a result of elastic interactions. Ex- 
amples include (2 x l)-0 domains on Cu(llO) 0], or- 
dered two-dimensional (2D) islands of Ag on Pt(lll) 0], 
stripes of alternating dimer direction on Boron doped 
Si(OOl) 0,0, square patterns of N/Cu(100) 0, stripe 
domains of Pb/Cu(lll) [l], and Au stripes on W(llO) Q. 
In all these cases, the length scales of the periodic modu- 
lations range from a few to hundreds of nanometers, well 
above the atomic distances. This is understood within 
the analytical theory Q , which states that the period of 
a given stripe pattern depends exponentially on the ratio 
of the short-range (boundary) and long-range (elastic) 
interaction energies. The period D can be written as 



D = 2iTa 
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where a is a microscopic length, Ci is the formation en- 
ergy of the stripe boundary, and C2 is the prcfactor of 
the energy gained by the elastic relaxation due to the 
formation of stripes The prcfactor of the exponent 
comes from a smearing of the stripe boundaries, and so 
is related to the boundary width. This result is identical 
to that obtained for magnetic layers with dipolar interac- 
tions 13, U I • We note that the correspondence between 
the dipolar Ising lattice and the 2D lattice gas is one- 
to-one, as the elastic interaction energy between lattice 
defects scales as 1/r^ [l^ . 

Studies of mesoscopic patterns show that the period- 
icity varies strongly with temperature. Moreover, there 
exists a transition temperature, at which the stripe phase 
disappears. This temperature dependence has been at- 
tributed to the reduction of the boundary free energy 
through density fluctuations as the disordering tempera- 
ture is approached [l3| • Earlier work on Ising spin lattices 



showed that the effect of fluctuations can be introduced 
by an entropic term, which results in a temperature de- 
pendent period for the magnetic stripes 14 1. 

The density fluctuations, which reduce the boundary 
free energy, also increase the boun dary (or domain wall) 
width with increasing temperature [l5| . Gehring and Ke- 
skin took this broadening into consideration, and they 
argued that the stripe disordering transition takes place 
when the boundary width becomes half the stripe pe- 
riod 0. 

All studies up to date treat the temperature regime 
close to the transition separately, pointing to a 
"crossover" as the boundary width becomes cornparable 
to the stripe period (a discussion can be found in [lO[). In 
this paper, we report the observation of a stripe phase in 
sub-monolayer Pd films grown on W(llO). We demon- 
strate that the temperature dependence can be repro- 
duced throughout the whole temperature range by a scal- 
ing of the boundary energy and boundary width. Us- 
ing low-energy electron microscopy (LEEM), we show 
that the triangular step-flow patterns observed be- 
low 1100 K turn into stripes that run along the [110] 
direction above this temperature. The period of the pat- 
tern sharply decreases from 200 nm to 57 nm with in- 
creasing temperature until the stripes disorder at around 
1170 K. The dependence of the period on Pd coverage is 
much weaker contrary to the expectation from the theory 
with sharp boundaries (a similar behaviour was found for 
Au/W(110) 0). 

The thermodynamics of Pd/W(110) can be described 
as that of a lattice gas with attractive interactions, and 
shows a condensate-gas transition with increasing tem- 
perature [3, [1^. Structurally, a Pd monolayer on 
W(llG) shows one-dimensional pseudomorphism along 
the [001] direction, whereas along [110] extra low-energy 
electron diffraction (LEED) spots are observable even at 
the highest temperatures corresponding to the Pd lat- 
tice relaxed through periodic dislocation lines [I^. This 
indicates a stress anisotropy with a low surface stress 
along [110], readily explaining the direction of the ob- 
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FIG. 1: (Color online) Monolayer Pd stripes on W(llO) with increasing temperature. The stripe direction is along [lIO]. The 
condensed Pd islands appear bright at the electron energy used (5.2 eV). The image size is 1 /im^. A 350 nm profile across the 
stripes is shown below each image. Note that the rectangular profile of the sharp boundaries is broadened due to the resolution 
of the microscope. In all images, monoatomic substrate steps with circular shapes can be recognized. 



served stripes in our study. 

The LEEM measurements were performed with the 
SPELEEM microscope at Elettra (Italy) The in- 

strument allows real-time monitoring of the growth and 
evolution of adsorbate layers with a lateral resolution of 
12 nm. The W(llO) sample was cleaned by annealing at 
1000 C in 2 X 10~^ mbar oxygen, followed by repeated 
high temperature flashes in ultrahigh vacuum to remove 
oxygen. The absence of oxygen was confirmed by the 
sharp, low background (1x1) LEED pattern. Pd was 
deposited on the tungsten surface using an e-beam evap- 
orator at a rate of 0.2 monolayer per minute. The rate 
was calibrated via the time needed to complete the first 
monolayer at 800 C. 

A series of LEEM images of alternating Pd stripes in 
condensed and lattice gas phases are shown in Fig. [T] 
The period of the pattern and the contrast between the 
two phases sharply decrease with increasing temperature 
up to about 1170 K, at which point the stripes fully dis- 
order. The variation of the periodicity is facilitated by 
the high mobility of the condensed regions. This leads to 
stripe fluctuations around the reversible stripe disorder- 
ing transition. 

The natural logarithm of the stripe period, D, is dis- 
played in Fig. [2] as a function of temperature. Ex- 
cept for temperatures very close to the stripe disorder- 
ing transition, the logarithm of the period shows a lin- 
ear trend, suggesting that the temperature dependence 
comes mainly from the energetic parameters in the ex- 
ponent in eq. ([T]). As we will discuss shortly, due to 
smearing of the stripe boundary, the period levels off at 
a nonzero value close to the disordering transition. 

As mentioned above, studies on the dipolar Ising lat- 
tice arc informative also for the behaviour of adsorbate 
stripes [2lj. Recent studies show that very close to the 
disordering transition, the stripe profile can be approx- 
imated by sinusoidals (keeping the lowest orders in a 
Fourier expansion) [2^ . The result is a quadratic de- 



pendence of the period on temperature obtained by min- 
imizing the free energy: 
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where Z?o is the period at the disordering temperature, c 
is a constant given as a combination of the interaction pa- 
rameters, and Tc is the disordering temperature. Despite 
being a good approximation close to the disordering tem- 
perature, it breaks down away from the transition due 
to the implicit assumption that the width of the stripe 
boundary is comparable to its period. Indeed, Fig. [1] 
shows that at the lower temperatures the stripe period 
is much larger than the boundary width, and that the 
two phases (i.e. dark and bright stripes) have unequal 
widths. Thus the lowest order components of a Fourier 
expansion are no longer sufficient to describe the stripe 
profile. 

Instead, the temperature dependence can be obtained 
from general scaling laws governing the parameters in 
eq. 11]). For the dipolar Ising lattice, the changes in 
the boundary free energy and boundary width have been 
identified as the determining factors |16| . The boundary 
free energy was found to decay linearly with temperature 
due to the increase in entropy [3, [Is] . The result is a 
scaling of the exponent in the stripe period: 
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where Ci is the temperature dependent boundary free 
energy, and is the transition temperature for the Ising 
model without the long-range dipolar or clastic interac- 
tions. On the other hand, the boundary width scales as 
the correlation length [16[. Using the critical exponent 
J/ = 1 for the 2D Ising model. 
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the transition temperature is given by the ratio of the 
long-range and short-range energies: 



FIG. 2: Natural logarithm of the stripe period as a function 
of temperature. The data points are shown as solid circles. 
The line corresponds to the fitting function in eq. ((Sjl. 



where we define the unsealed boundary width as wq = 
2TTa. Combining eqs. ([T]), ([3]) and we can write the 
logarithm of the stripe period as 



ln{D) ^ 
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The basic features in Fig. [5] can now be understood. The 
first term on the right hand side of eq. ([5]) dominates at 
the lower temperatures giving the linear trend in ln{D) 
vs T . The second term has a mild temperature depen- 
dence away from the transition, but becomes increasingly 
important close to it. This corresponds to a broadening 
of the stripe boundary, which slows down the variation 
in the period and levels it off at the stripe disordering 
temperature. 

The resulting fit to the experimentally determined 
stripe period is shown as a solid line in Fig. [21 The fit- 
ting parameters were found to be wq — 1.95 ± 0.04 A, 
r° = 1201.7 ± 2.6 K and the energy ratio C1/C2 = 
39. 8 ±0.8. The scaling gives excellent agreement with the 
data throughout the whole temperature range for which 
stripes are observed. It should also be noted that near the 
disordering temperature, Tc, the quadratic dependence in 
eq. ([2]) can be retrieved from the functional form used for 
the fit (see Appendix). 

The value of T° reveals an important effect. The scal- 
ing laws in eqs. ^ and ^ were found from the 2d Ising 
model. They do not take into account the long-range 
elastic interactions, which counteract the short-range in- 
teractions between Pd atoms and favor disorder, lower- 
ing the transition temperature from the pure Ising value, 
T°, to the stripe disordering temperature, Tc = 1170 K. 
Moreover, as predicted 2J], we see that the decrease of 
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In physical terms, the transition temperature is lower 
than the value expected from a diverging correlation 
length, because the disorder sets in when the boundary 
roughness (or the correlation length) becomes compara- 
ble to the stripe period. The same phenomenon is ob- 
served in finite size Ising lattices, in which the disorder 
transition takes place when the fluctuations reach the 
system size. Such a reduction in the transition tempera- 
ture was shown to scale inversely with the finite size [2^ . 
The period of Pd stripes at the transition is approxi- 
mately given by Dq = e^wo(Ci/C2) from Eqs. (O and 
([5]), providing an analogous relation between the period 
and the lowering of the transition. 

The actual condensate to lattice gas phase transition, 
which should take place at r°, can be observed in the 
intensity of the Pd diffraction spots shown in the inset 
of Fig. [21 In agreement with the discussion above, we 
observe that the crystalline order of the Pd layer persists 
after the stripe contrast vanishes at Tc. A fit to the 
temperature dependence of the diffraction intensity (see 
Fig. [3]) gives a condensate to lattice gas transition at 
T° = 1210 ± 5 K, only slightly higher than the value 
extracted from Eq. ([S]). The intermediate phase, between 
Tc and Tg , does not have the long-range stripe order, 
however it consists of condensed islands buried in the 



lattice gas phase as discussed in [25 1. 

The exponential decay of the Pd diffraction spot inten- 
sity with temperature up to about 1170 K is much faster 
than that expected for a physically realistic Debye- Waller 
factor. It is rather due to the decrease of the number of 
atoms in the crystalline phase by sublimation into the 
lattice gas. Consequently, the energy in the exponen- 
tial, Ek = 0.88 eV, corresponds to the removal of Pd 
atoms from the kink sites of the stripe boundaries. This 
translates into a boundary energy of Ci = 197 meV/A 
along [110]. The value of the boundary energy, along 
with C1/C2 = 39.8, results in an elastic energy parame- 
ter C2 ~ 5.0 meV/A. Using the bulk elastic constants for 
tungsten, we obtain the surface stress change across the 
Pd step to be Arl^"^' « 4.6 N/m (assuming an isotropic 
elastic response for the tungsten substrate according 
to [2^). This is in very good agreement with density- 
functional theory calculations, which suggest that the 
tungsten surface stress along [001], 5.26 N/m 27j . is re- 
duced nearly to zero upon Pd adsorption [28| . 

Finally, we note that no orientational melting of the 
stripes is observed. With increasing temperature, con- 
trast between the condensate and gas stripes decreases 
monotonically, until it disappears at Tc. However, the 
direction of the pattern is preserved all the way up to Tc 
as shown in Fig. [1] 
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where Dq is the period at Tc, and T — T^ — ST. 
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FIG. 3: The intensity of the Pd extra LEED spot (shown by 
the arrow in the inset) as a function of temperature. Electron 
energy is 35 eV. The solid line is a fitting function of the form 
gBfc/fcTj--j^ _ T/T^")^'', where 13 = 1/8 is the critical exponent 
for the order parameter in the 2D Ising model. Ek is the 
energy of sublimation from the condensate to the lattice gas 
phase. 



In conclusion, we have shown that stress-induced 
stripes of Pd on W(llO) behave according to the scaHng 
laws for the 2D Ising model with increasing temperature. 
Using low-energy electron diffraction and microscopy, we 
have demonstrated that the condensate to lattice gas 
transition takes place at a slightly higher temperature 
than the disordering of the stripes. In analogy with the 
lowering of the transition temperature in the finite-size 
Ising lattice, this difference in the disordering temper- 
ature of the atomic and mesoscopic scales is explained 
by the loss of stripe order when the correlation length 
becomes comparable to the stripe period. 

We thank Natasa Stojic for making available the re- 
sults of her surface stress calculation on Pd/W(110) be- 
fore publication. 

Appendix. Eq. ^ can be obtained by expanding 
eq. ([5]) near the stripe disordering temperature, Tc. We 
first rewrite the scaling expression for a temperature 
slightly below Tc. 



D{Tc - ST) = 
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Using eq. Q we can expand both the denominator and 
the exponent in powers of (Ci/C2)((5T/T°). The result 
has the same form as in eq. ([5]): 
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